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Abstract 

We compute the entropy of entanglement in the ground states of a general class 



■ of quantum spin-chain Hamiltonians — those that are related to quadratic forms of 
. Fermi operators — between the first N spins and the rest of the system in the limit 
I of infinite total chain length. We show that the entropy can be expressed in terms of 
' averages over the classical compact groups and establish an explicit correspondence 

between the symmetries of a given Hamiltonian and those characterizing the Haar 
C^! measure of the associated group. These averages are either Toeplitz determinants 

' or determinants of combinations of Toeplitz and Hankel matrices. Recent gener- 

■ alizations of the Fisher-Hartwig conjecture are used to compute the leading order 
^ , asymptotics of the entropy as ^ oo. This is shown to grow logarithmically with 

N. The constant of proportionality is determined explicitly, as is the next (constant) 
term in the asymptotic expansion. The logarithmic growth of the entropy was previ- 
^ , ously predicted on the basis of numerical computations and conformal-field-theoretic 

I calculations. In these calculations the constant of proportionality was determined in 

terms of the central charge of the Virasoro algebra. Our results therefore lead to an 
explicit formula for this charge. We also show that the entropy is related to solutions 
of ordinary differential equations of Painleve type. In some cases these solutions can 
be evaluated to all orders using recurrence relations. 
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1 Introduction 

Random Matrix Theory, as developed by, amongst others. Freeman Dyson, provides the 
natural framework for calculating statistical properties of quantum fluctuations. It has 
had a profound impact in many of the traditional application areas of quantum mechanics, 
including condensed matter physics and nuclear physics. Our purpose here is to establish 
a new application: to calculating entanglement in quantum spin chains. 

The importance of entangled quantum states lies in their ability to exhibit correlations 
that cannot be accounted for classically. This feature of quantum mechanics has been 
known to physicists for almost seven decades — since Einstein, Podolsky and Rosen in- 
troduced their famous gedankenexperiment. However, it is only relatively recently that 
entanglement has come to be viewed as a physical resource for manipulating quantum in- 
formation. As for any other physical resource, like energy or entropy, it important to be 
able to quantify entanglement; that is, to assign a measure to it. This aspect of entangled 
states is still poorly understood, especially when the entanglement is shared between more 
than two systems. However, when the entanglement of a pure state is shared between two 
parties, i.e. in a bipartite system, Bennett et al P have shown that it is consistent to define 
it as the von Neumann entropy of either of the two parts. 

We consider here the general class of quantum spin chains arising from quadratic chains 
of fermionic operators in their ground state. These systems are partitioned into two con- 
tiguous subchains. If the ground state is non-degenerate, this subdivision creates a pure 
bipartite system and we investigate its entanglement of formation. 

As is well known, the systems we are studying exhibit quantum phase transitions. 
These manifest themselves as qualitative changes in the decay of correlations: algebraic 
in the proximity of a critical point and exponential decay away from it. Entanglement 
plays a fundamental role in the quantum phase transitions that occur in interacting lattice 
systems at zero temperature |21 El IH E] • Under these conditions the system is in the 
ground state, which is also a pure state, and any correlations must be a consequence of the 
fact the ground state is entangled. It follows immediately that the entanglement changes 
qualitatively in the proximity of critical points. 

Vidal et al jS] studied the ground states of the XY and XXZ models, partitioning them 
into two consecutive subchains; they observed numerically that, when the Hamiltonian 
undergoes a phase transition, the entanglement of formation of this bipartite system grows 
logarithmically with the size N of one of the two parts. Jin and Korepin [3] considered 
the XX model, which corresponds to an XY spin chain with an isotropic interaction, 
and expressed the von Neumann entropy in terms of a Toeplitz determinant, which they 
evaluated asymptotically using the Fisher-Hartwig conjecture [Zj. Their result coincides 
with the numerical observations of Vidal a/ in that they found that the entropy 
grows like | logg A^. Subsequently Korepin ^ and Calabrese and Cardy P showed, using 
conformal field-theoretic arguments, that the logarithmic divergence of the entanglement in 
one dimensional systems is a general consequence of the logarithmic growth of the entropy 
with the size of the system at phase transitions. These arguments determine the constant 
multiplying the leading order logg A^ term in the asymptotics to be one-third of the central 
charge of the associated Virasoro algebra. 

Our goal here is to show that if a quantum spin-chain Hamiltonian posses certain 
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symmetries, the entanglement can be expressed as an average over one of the classical 
compact groups equipped with Haar measure, i.e. one of the following groups: U(A^), 
Sp(2A^) and O^(A^), where the superscript ± indicates the connected component of the 
orthogonal group with determinant ±1. Indeed, we establish a precise correspondence 
between the functional form of the appropriate Haar measure and the symmetries of the 
Hamiltonian. The XX model turns out to be an example of a system with U(A^) symmetry. 
The averages that occur can be expressed either as Toeplitz determinants, in the case of 
U(iV), or as determinants of specific combinations of Toeplitz and Hankel matrices for the 
other compact groups. Recently, Basor and Ehrhardt |H] and Forrester and Frankel [H] have 
computed asymptotic formulae for such determinants, and these allow us to write down 
the leading-order and next-to-leading-order terms in the asymptotics of the entanglement 
in the limit as the total number of spins tends to infinity and then as ^ oo. 

We find that in the proximity of a critical point the entanglement grows logarithmically 
with A^, in agreement with the prediction of Korepin jS] and Calabrese and Cardy |E] . We 
derive a general formula for the associated constant of proportionality. This is a rational 
number, the numerator of which is shown to factorize into a universal part, related to sym- 
metries of the quantum Hamiltonian and which can be calculated from the random-matrix 
averages, and a non-universal (i.e. Hamiltonian-specific) part, which we also evaluate. Com- 
paring with Korepin's result leads to an explicit formula for the central charge, which plays 
a fundamental role in the conformal-field-theoretic approach. In addition, we obtain the 
sub-leading (constant) term in the asymptotics. 

We also show that the random-matrix averages are related to solutions of certain 
Painleve equations. In the case of Hamiltonians possessing U{N) symmetry, these so- 
lutions, and thus the entropy of the entanglement, can be determined exactly by means of 
recurrence relations. 

The structure of the paper is as follows. In section |21 we review the results of Vidal et 
al [3] and Jin and Korepin |@]. In section El we introduce the models that we investigate. 
Sections IH El El and [3 are devoted to identifying the symmetry classes of Hamiltonians 
that are associated to each compact group. In section El we apply generalizations of the 
Fisher- Hart wig conjecture to compute the entanglement for each symmetry class. Finally, 
in section El we investigate the connection between the interpretation of the entanglement 
as an average over a compact group and the theory of Painleve equations. 



2 Bipartite entanglement in the XY model 

We begin by reviewing the results of Vidal et al |3j and Jin and Korepin |lj concerning the 
entanglement in the ground states of the XY and XX models. These systems are spin-1/2 
ferromagnetic chains with an exchange coupling a in a constant transversal magnetic field 
h. The Hamiltonian H = hHa with given by 

A/-1 M-l 

Ha = ~Y.[(^+ ^^l^U + (1 - ^y>U - E (2-1) 

where a"" denotes the Pauli matrices and a = x,y, z. It is sometimes convenient to work 
with the spin operators Sj = o'j/'2 instead of the Pauli matrices. The parameter 7 lies in 
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the interval [0, 1] and measures the geometric anisotropy of H^- Throughout this paper we 
will assume periodic boundary conditions, i.e. a%j = ctq. At zero temperature the system is 
in the ground state |^'g). In the limit M — >■ oo the Hamiltonian ()2.H1 undergoes a quantum 
phase transition at = l/{2h): the spin-spin correlation functions {^g\SjS^ |^'g) and 
(^'gl SjSl l^'g) decay algebraically at the critical point and exponentially when a ^ a^- 

Let us partition the Hamiltonian ()2.1|) into two subchains, which we denote by P and 
Q. The subsystem P is composed of the first iV spins, with 1 -C -C M. The Hilbert 
space where the Hamiltonian ()2.1|1 acts is the direct product H. = Hp ^Hq, where Hp and 
Hq are generated by the spins in P and Q respectively and are spanned by the vectors 

N-l M-1 

Hia-y^l^p) and l[{a;Y^ \^p) , r, = 0, 1, (2.2) 

j=0 j=N 

where = {a^ ±iaj)/2 and |^'f) denotes the ferromagnetic state with all spins up. Since 
the ground state is non-degenerate this subdivision creates a pure bipartite system. The 
measure of entanglement is defined as the von Neumann entropy of either subchain: 

Ep = Eq = - Tt pp log2 Pp = - Tr pq log2 pq, (2.3) 

where pp and pq are the reduced density matrices of P and Q, i.e. 

Pp = Trqppq and pq = Trpppq. (2.4) 

The operators ppq form the density matrix of the whole system, ppq = |'I'g)(^'g|. 

Let us introduced the Jordan- Wigner transformations at each site of the lattice {1, . . . , M}: 

m2z+i = ^ JJ j and ma; = ^JJajj af. (2.5a) 



The inverse relations are 



ai = 11712117121+1, 

erf = I JJim2jm2j+i I 17121+1, 

\j=0 

o-f = I JJ im2jm2j+i j 17121 



(2.5b) 



These operators are Hermitian and obey the anticommutation relations {mj,7nk} = 26jk- 
We also define 

bi = {m2i+i - im2i)/2 and b] = (m2/+i + im2i)/2, (2.6) 
which obey the anticommutation relations 

{b,,bk} = and {b],bk} = Sjk (2.7) 
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and therefore are Fermi operators. 

The expectation values of the operators ()2.5a|l with respect to l^'g) are 

(*g|mfc|*g) = 0, (2.8a) 
(*g| mjuik |*g) = 5jk + KCM)jk, (2.8b) 

where the correlation matrix Cm has the block structure 



C, 



M 



with 



For large M, the real numbers gi are the Fourier coefficients 



/Cn 


C12 


CiM^ 


C21 


C22 


C2M 


\Cmi 


Cm2 


■ ■ ■ Cmm J 









jk = 


-9k-j 


)- 



(2.9a) 



(2.9b) 



g, = ^ r acosg-l + ia7sing ^^.lO) 
^"^■^0 ^/(ttcos^- l)^ + 72a2sin2^ 



Formulae ()2.8p were first computed (in a slightly different form) by Lieb et al [TU] and 
Barouch and McCoy jTI]. The expectation values of products of arbitrary numbers of the 
operators ()2.5a|l can be obtained from ()2.8j) using Wick's theorem. Equation ()2.8a|l follows 
from the invariance of under the map (aj, aj, cr|) y-^ ""'^j ' "^i)' appendix B we 

extend ()2.8b|l to a class of Hamiltonians that generalizes the XY model. 

In order to compute the entropy ()2.3|1 we need an expression for the reduced density 
matrix pp. We therefore restrict our analysis to the first consecutive spins. The correla- 
tion matrix is simply constructed by removing the last 2(M — A^) rows and columns from 
the matrix ()2.8b|) : we shall denote this restriction by C^. If V is an orthogonal matrix, 
then the operators 

2Af-l 

dj = ^ Vjkmk (2-11) 

fc=0 

are Hermitian and obey the anti-commutation relations {dj,dk} = 26jk too. Furthermore, 
since the Fourier coefficients defined in equation ()2.10|) are real, there exists an appropriate 
V G S0(2A^) that block-diagonalizes C^'- 

VC^V'J^uJ^^ J), (2.12) 

where the uj are real and, for reasons that will become apparent later, lie in the interval 
[— 1, 1]. Now, let us define 



Cj = {d2j+i - id2j)/2, j = 0,...,N -1, 



(2.13) 
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which are Fermi operators, i.e. 

{cj,Ck} = and {c],Ck} = 5jk (2.14) 

and are hnearly related to those ones defined in ()2.(jj) . A basis for Tip is also given by the 
2^ vectors 

N-l 

n(c5)'-M*vac), r, =0,1, (2.15) 

j=0 

where the state l^'vac) is defined by the condition 

C,|*.ac)=0, j=0,...,iV-l. (2.16) 

The expectation values of the Cjs are readily obtained from equation ()2.12j) : 

(*g|c,|*g) = (*g|c,Cfc|*g) = 0, (2.17a) 
(*j4c,|*g) = 5,fci^. (2.17b) 

The reduced density matrix pp can be computed directly from the expectation val- 
ues (I2.17j] : we describe the details of the derivation in appendix A. The final expression is 

P^-li(^4-' + ^^'4)- (2.18) 

i=o ^ ^ 

In other words, as equations ()2.17|) already suggest, these fermionic modes are in a product 
of uncorrected states, therefore the density matrix is the direct product 

N-l ^ _ ^ ^ _^ ^ 

pp = (g)p, With p^ = ^^c]c, + ^^c,c]. (2.19) 

j=0 

As a consequence, if one of the ujs lies outside the interval [—1, 1], then either (1 + z/j)/2 
or (1 — z/j)/2 would be negative and so pj could not be a density matrix. At this point 
the entropy of the entanglement between the two subsystems can be easily derived from 
equation ()2.3|1 : 

Af-l 

Ep = 5^e(l,z/,), (2.20) 

j=0 

where 

X + ly ^ f X + h'\ X — V ^ f X — v\ 
e(., .) = -—log, (_ j - _ log, (— j . (2.21) 

Vidal et al 0] computed Ep numerically for the XY model and observed that for values 
of a and 7 close to a critical point it grows logarithmically with A^, while away from the 
phase transitions it is either zero or saturates to a constant value. Jin and Korepin |3] 
computed Ep when 7 = 0, i.e. for the XX model, using the Fisher- Hartwig conjecture, 
which gives the leading order asymptotics of determinants of Toeplitz matrices whose 
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symbols have zeros or discontinuities. This is also central to our own approach, so we will 
give a brief account of their method. 

When 7 = the numbers ()2.1()|1 are the Fourier coefficients of the step function 

jl if-k<e<k 
^ ~ [-1 iik <e <2ti - k, ^ ' ' 

where k = arccos(l/Q;). In order to be in a critical regime, the parameter 1/a must lie in 
the interval (—1, 1); for |a| < 1 the entanglement Ep is trivially zero. In physical terms this 
means that outside the critical regime the magnetic field is strong enough to align all the 
spins, therefore entangled states cannot appear. Since g{6) is even its Fourier coefficients 
have the symmetry gi = g^i, therefore the correlation matrix Cn factorizes into the direct 
product 

CN = T^[g]®(^\ J), (2.23) 

where T]sf[g] is the matrix 

{TN[g]),k = gj-k. (2.24) 

Matrices of this type are called Toeplitz matrices and g{6) is called the symbol of Tisi[g]. 
Toeplitz matrices and their determinants, known as Toeplitz determinants, will play an 
important role in our analysis. 

As a consequence of ()2.23j) the UjS are simply the eigenvalues of T^lg]. One can then 
rewrite ()2.2()j) using Cauchy's theorem j^: 

Ep = lim lim ^ / e(l + e, X)^^^^^0^dX, (2.25) 

where Dj\f{X) is the characteristic polynomial 

DN[g]{X) =det {XI -T^lg]), (2.26) 
which is also a Toeplitz determinant with symbol 

9{0-A) = S ~ 2.27 

^ ^ |A + 1 if A;<^<27r-A;. ^ ^ 

The contour of integration c(e, 6) depends on the parameters e and 6 and includes the 
interval [—1, 1]; as e and 6 tend to zero the contour approaches the interval [—1, 1]. This 
technical expedient is necessary in order to guarantee that the branch points of e(l + e. A) 
lie outside the contour of integration and thus that e(l + e. A) is analytic inside c(e, 6). As 
already mentioned, the polynomial Dn{X) can be estimated for large N using a formula, 
which we will discuss in detail in section known as the Fisher-Hartwig conjecture [Zj . 
This asymptotic formula can then be inserted in equation ()2.25|1 and the integral explicitly 
evaluated. We will not go into the details of the computation, which can be found in Jin 
and Korepin's original paper. Their result is that for the XX model 

Ep log2 A^, A^ ^ oo. (2.28) 
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This logarithmic divergence is a direct consequence of a general theorem on the spectra 
of Toeplitz matrices jT2I, which in our case reduces to the statement that in the limit 
oo all except for O(lniV) of the eigenvalues of T^lg] migrate toward the points 1 and 
— 1, where, trivially, we have e(l, 1) = e(l, —1) = 0. 

We now make an observation relating to ()2.25|) that will have important consequences 
for the way we shall proceed. Toeplitz determinants are equivalent to averages of appropri- 
ate functions over the group of unitary matrices U(iV). Let us make this statement more 
precise. The group U{N) is compact and therefore has a unique left and right invariant 
measure known as Haar measure. The invariance of the measure makes it a natural proba- 
bility density for unitary matrices, because different regions in U(A^) are equally weighted. 
This is usually referred to as the Circular Unitary Ensemble (CUE) of random matrices. 
After having integrated over the degrees of freedom associated with the eigenvectors, the 
explicit expression of Haar measure becomes 



Pt 



V{N) 



{6i, . . . , On) 



1 



(27r)^A^! 



n 



(2.29) 



^<j<k<N 



Now, let G{U) be a function on U(A^) that depends only on the eigenvalues of U and is 
such that 



N 



G{U) = \{g{ 



(2.30) 



where g{9) is 27r-periodic. An identity due to Heine ^3] and — in the form that we use 
Szego [H], pp. 27 and 288, asserts that 



G{U) 



\]{N) 




(2vr)^Ar! J, 

det (fi'i-fc)j_jt=o,...,Ar-l ' 



n 



d^l ■ ■ ■ d^ 



(2.31) 



N 



^<j<k<N 



where 



9i 



^ / 9(9) e-^^' 
2vr Jo 



de. 



(2.32) 



Throughout this paper the brackets (■) denote the average with respect to Haar measure 
(not necessarily just that of U(A^)). 

For the XX model this remarkable identity allows us to express ()2.25j] as an average over 
U(A^). Two questions now arise. First, can ()2.25j) be generalized to other Hamiltonians? 
And second, since there are no obvious reasons why U(A^) should be singled out with respect 
to the other compact groups, are there systems for which the symmetries of the interaction 
give rise to a von Neumann entropy that can be re-expressed in terms of averages over 
Sp{2N) and O^(A^)? 
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We note first that when 7 7^ the symbol of the matrix Tfyf[g] is 

a cos 9 — 1 + ia'-y sin 9 



'a cos 9 — Vf + ^"^a^ sin^ 9 



(2.33) 



The numbers Uj that appear in formula ()2.12|1 are not the eigenvalues of Tj^[g\, which is 
not symmetric and as a consequence need not have a real spectrum. Instead they are 
the eigenvalues of the matrix S = {T]\f[g]TN[gY)^^'^ . As a consequence, the entanglement 
is not straightforwardly expressible in terms of a Toeplitz determinant and so in terms 
of an average over U(A^). Even so, it can again be expressed as an integral transform 
of the characteristic polynomial of S. It is worth remarking that the spin-spin correlation 
functions can still be expressed in terms of Toeplitz determinants and so related to averages 
over U(A^). 

In the following sections we shall show that for Hamiltonians that can be expressed 
as quadratic forms of fermionic oscillators, of which the XY model is an example, there 
exists a one-to-one correspondence between the symmetries of the Hamiltonian and the 
interpretation of the von Neumannn entropy in terms of averages over one the classical 
compact groups when 7 = 0. In sectionlHlwe establish a connection between a generalization 
of formula ()2.25|1 and the theory of Painleve equations. 



3 Quadratic chains of fermionic operators 

The most general form of Hamiltonian related to quantum spin chains is 



a 



A/-1 



2 y: h% 



(3.1: 



j=0 



where a and 7 are real parameters, < 7 < 1 and the fo^s are the Fermi oscillators defined 
in equation ()2.6|) . We take periodic boundary conditions, i.e. 6^ = h^. Since is 
Hermitian, A must be a Hermitian matrix, and because of the anticommutation relation of 
the Fermi operators, B must be antisymmetric. Without loss of generality, we will consider 
only matrices A and B with real elements. 

Up to a overall constant, the XY model 1)2.111 maps into the Hamiltonian 1)3.1)) '^ with 
the matrices A and B given by 




and B 



/O 







1 
1 
-1 















-1\ 
























1 

0/ 



(3.2) 



^This is strictly true only for open-end Hamiltonians. If we assume periodic boundary conditions, then 

bM-ibo- However, because we are interested 



the term b\j_^bo should be replaced by Il^io^ ^ ^ 

in the limit M 00, the extra factor in front of 6a/-i^o can be neglected. 
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Similarly, equation ()2.5b|) maps the Hamiltonian ()3.H) into the spin chain 

/ k-l 



0<i<fc<A/-l 



izr ' (3.3) 



In other words, the two systems ()3.Hl and ()3.3|1 are equivalent. One sees that the second 
sum in ()3.H) is the analogue of a uniform magnetic field, the parameter 7 introduces a 
geometric anisotropy in the interaction, and a is an exchange coupling constant. 

We will here be concerned with the entanglement between the first N oscillators and 
the rest of the chain, when the system is in the ground state |^'g) and as the length of chain 
tends to infinity. In a similar fashion as for the XY model we decompose the Hilbert space 
into the direct product Ti = Hp ® Hq, where Hp is generated by the first sequential 
oscillators and Hq by the remaining M — N. They are spanned by the vectors 

N-l M-N 

n(&5PI*vac) and l[{b]Y^\^..c), = 0,1, (3.4) 

j=0 j=N 

respectively, where the vacuum state l^'vac) is defined by 

6,|*vac)=0, j = 0,...,M-l. (3.5) 

Our goal is to determine the asymptotic behaviour for large A^, with = o(M), of the 
von Neumann entropy 

Ep = - Ti pp log2 Pp, (3.6) 

where pp = Trqppq and ppq = |*g)(*g|. 

The first step involves determining the expectation values with respect to |^'g) of prod- 
ucts of arbitrary numbers of operators (j23aj). From the invariance of the Hamiltonian ()3.1|) 
under the transformation bj ^ —bj, it follows that {^g\mi |^'g) = 0; for the same rea- 
son, the expectation value of the product of an odd number of rrijS must be zero. The 
expectation values of the product of an even number of rrijS can be computed using Wick's 
theorem 

g\ mjj^mj^ ■ ■ ■ TUj^^ \^ g) = {~^y Y\. (contraction of the pair) , (3.7) 

all pairings all pairs 

where a contraction of a pair is defined by (^'g| rrij^mj^ I'^g) and p is the signature of the 
permutation, for a given pairing, necessary to bring operators of the same pair next to one 
other from the original order. 

Before continuing our analysis it is worth noting that the spin-spin correlation functions 
(*g| S^S^ |*g) and (*g| SljSl |*g) of the XY models are products of the type These 
correlations were first studied by Barouch and McCoy ^T], who showed that in general 
they are Pfafiians, which, when the system reaches thermal equilibrium, reduce to Toeplitz 
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determinants and can be computed using Szego's theorem and the Fisher-Hartwig conjec- 
ture. The different behaviour of these determinants away from and in the proximity of 
critical points determines how quantum phase transitions appear. Thus, in the same way 
that we can ask which features of the Hamiltonian p.ip lead to expressions for the entropy 
of entanglement in terms of an average over \J{N), so we can enquire which properties of 
Ha ensure that its physical correlations are expressible as averages over U(iV). In view of 
the way physical correlations and entanglement affect each other, these two questions are 
closely related. Similarly, we can pose the same questions for the other classical compact 
groups Sp{2N) and 0^{N). 

The expectation values (^'g| mjmk |^g) can be deduced from the work of Lieb et al |l()j : 
as already mentioned in section |2l we give the details of this computation in appendix B. 
They generalize formula ()2.8b|) . in that the Fourier coefficients Qj^k and Qk-j in the 2x2 
block ()2.9b|) are replaced by the matrix elements {TM)jk and (TM)kj, where the matrix Tm 
is defined by 

M-l 

(^M),fc = 5Z ^^.-0^^, 3,k = 0,...,M-l, (3.8) 

1=0 

and the vectors 0^ and -j/'a; are real and orthogonal and obey the eigenvalue equations 

a^{A - {2/a)I - -iB){A - i2/a)I + 7^)0^ = |Afc| Vfc, (3.9a) 
a\A - {2/a)I + ^B)iA - {2/a)I - ^B)^/,, = |Afc| (3.9b) 

These vectors are related by 

a(A-(2/«)/ + 75)<^, = lAfclV^fc, (3.10a) 
aiA - {2/a)I - ^B)iP, = |Afc|0,. (3.10b) 

From equations ()3.9j) it follows that the 0^s and V'fcS are eigenvectors of positive symmetric 
matrices, therefore the corresponding eigenvalues are always real and positive; for later 
convenience, here and in appendix B we express them as the square of the absolute value 
of a complex number. 

Now, the derivation of the von Neumann entropy ()3.6|) is identical to the analogous 
computation for the XY model described in section |21 The formula for the entropy of the 
subchain P that one obtains is 

Ep= hm hm — / e(l + e,A)^^^^^^dA, (3.11) 

where 

Dn{X) = det {IX- S), (3.12) 

S is the real symmetric matrix (T/vT^)^^^ and T/v is obtained from the matrix ()3.8j) by 
removing the last M — N rows and columns. The path of integration of the integral 
in p.lip is the same as that in ()2.25|) . In order for p.llj) to have physical meaning, the 
eigenvalues of S must all be in the interval [—1, 1]. This will have to be verified for the 
various Hamiltonians of the type ()3.1|1 that we shall consider. 
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4 U(A^) symmetry 

If there is a connection between the entanglement of the ground state or, more generally, 
between the physical correlations of the ground state and the classical compact groups, it 
must be reflected somehow in the symmetries of the system of fermionic oscillators (jH.lj) . 
The most obvious symmetry that the XY model has in common with U(iV) is the invariance 
under translations: for the XY model it manifests itself in the fact that the matrices ()3.2|) 
are mapped into themselves if their rows and columns are simultaneously shifted by the 
same integer q (the periodicity of Ha is of course inherited by the periodicity of the rows 
and columns of A and B); for U(A^) such symmetry is reflected in the invariance of the 
integral ()2.HH) if all the 6jS are shifted by the same amount x- It is worth noting that the 
integral ()2.31|) can also be written in the form 




2n 




det [Sn (9, - 9k\,^,_^ d^i ■ ■ ■ d^,v, {^X 



U(iV) 

where Sn{z) is the kernel of the Haar measure of U(A^) and is given by 

^ ^ 2tx sin(2/2) ^ ' 

The invariance under translations of the Haar measure is then evident in the kernel ()4.2|) . 
We shall now determine a general formula for the matrix ()3.8j) in the limit M — > cx) by 
assuming that the physical interaction in the Hamiltonian (jH.lj) is invariant under transla- 
tions; we shall also deduce under which conditions the characteristic polynomial appearing 
in the integrand of formula ()3.11|) can be interpreted as an average over U(A^). 

In order to simplify the algebra we shall denote A = aA — 21 and B = ajB. If Ha 
is invariant under translations of the lattice {0, 1, . . . , M — 1}, than the elements of the 
matrices A and B must depend only on the difference between the row and column indices, 
i.e. A and B must be Toeplitz matrices. In addition, because of the periodic boundary 
conditions, A and B must be cyclic. 

Now, let a and b be two real functions on Z/MZ, even and odd respectively. The matrix 
elements of A and B can be written as 

A^k = a{j-k) and Bjk = b{j-k), j, = 0, . . . , M - 1. (4.3) 

A complete set of orthonormal eigenvectors of cyclic matrices are the complex exponentials: 



j,/c = 0,...,M-l, (4.4) 



as can be easily verified by direct substitution. Now, the matrices A and B defined in 
equation ()4.3|1 commute. This becomes straightforward if we notice that {AB)jk is the 



Random Matrix Theory and Entanglement in Quantum Spin Chains 



13 



convolution of a and b evaluated at j — k: 

M-1 M-1 



M-1 ^ (4-5) 

= ^ a(/ - k)bU -l) = {BA)jk. 



As a consequence, the complex exponentials ()4.4p are a complete set of eigenvectors of 
A + B too. Therefore, we have 



M-1 



J2 Hp - + Kp - 0] (Pki = K<Pkp = lAfcl^fcp, (4.6) 



1=0 



with 'j/'fc = 4>k-^'k/\^k\- Because both 0^ and 'j/'a: normalized, A'^/|Afc| is a phase factor 
and we can set A'^, = A^. It is important to notice that since we have taken the 4>i^s and 
V'fcS to be common eigenvectors of A and B they are not in general real. The matrix ()3.8|1 
should therefore be replaced by 



M-1 



= 5Z ^'i'^"^' J,fc = 0,...,M-l, (4.7) 

1=0 

which a priori need not be real. However, the matrix ()4.7p is unitarily equivalent to ()3.8p : 
indeed, we shall prove that ()4.7|) is real too and as consequence there is an orthogonal 
matrix O such that OTmO^ = Tm- Therefore, we shall not distinguish between them. 

The eigenvalues of A + B can be determined simply by inserting the eigenvectors ()4.4|1 
into the left-hand-side of equation (j4.(j|l and using the parities of the functions a{j) and 
b{j). We have 

^ ^ f a(0) + 2 Efir'^^' cos kj + ib{j) sin kj) if M is odd 

^ \a(0) + (-iya(M/2) + 2^^!/-f"^ (a(j)cosA;j +i6(j)sinA;j) if M is even, ^ ' ^ 

where k does not denote an integer anymore but the wave number 

Ott/ 

k = j^, / = 0,...,M-1. (4.9) 

The matrix ()4.7|) now becomes 

277(1- l/M) 

(TM),, = i- ^e-^(-')Afc, (4.10) 

fc=o l^^'^l 



where Ak = 2tt /M. From ()4.10p it follows that taking the complex conjugate of the right- 
hand side of ()4.10j) is equivalent to replacing A; by 27r — in the sum, therefore Tm is real. 
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If now we focus on the system consisting of the first oscillators and let M — >■ cxo, we 
obtain 

where A (6*) is the periodic function 

oo 

A(^) = J2 Aj-e'^'^ (4.12) 

j=-oo 

with Aj = a{j)—b{j) if j > and Aj = a{j)+b{j) if j < 0. (We have also implicitly assumed 
that as j — > oo, a(j) and b{j) tend to zero sufficiently fast for the Fourier series (j4.12|l 
to converge.) Thus, T]sf[g] is a Toeplitz matrix with symbol g{6) = A(6')/|A(^)|. It is 
worth emphasizing that formula ()4.1H) has been obtained by assuming only the translation 
invariance of the Hamiltonian ()3.H1 and periodic boundary conditions. Finally, if we define 

Mj) = I' if^-^'-^dAT ^^^^^^ 
I otherwise, 

then we recover the XY model with the choice 

a{j) = aMj)+5,{-j)]-25oij) and bU) = -a^ Mj) - 5,i-j)] . (4.14) 

At this point a few remarks relating to ()4.11|) should be made. First, because the inter- 
action is invariant under translations, we can take any set of A^ consecutive oscillators and 
the corresponding matrix T^lg] will still be a Toeplitz matrix. Second, as a consequence 
of the Jordan- Wigner transformations (j2.5bjl and Wick's theorem ()3.7|1 the spin-spin cor- 
relation functions are Toeplitz detrminants, i.e. averages over U(A^). Finally, there are 
important implications for formula ()3.11|) . We mentioned that in order to have physical 
meaning the eigenvalues of {T]\j[g]T]\f[gYy^'^ must lie in the interval [—1, 1]. Since the sym- 
bol g{d) has absolute value one, a theorem on the spectrum of Toeplitz matrices ^2] — the 
same theorem mentioned in section El — states that the eigenvalues of T^lg] are all inside 
the unit circle and approach the image of g in the limit N ^ oo. It follows that all the 
eigenvalues of {TN[g]T]sf[gYy^'^ lie in the interval [—1, 1]. It remains to establish when ()3.11|) 
is an average over U(A^). The condition is that T/v[5'] should be symmetric, in which case 
the correlation matrix Cn factorizes into the direct product as in ()2.23|1 . A necessary and 
sufficient condition in order for T^lg] to be symmetric is that A{6) should be real and even, 
or equivalently 7 should be zero; in other words, the interaction in the Hamiltonian ()3.1|) 
must be isotropic. When 7 = the symbol g{9) is a piece- wise continuous function that 
takes the values 1 and —1 and has discontinuities at all points 6r where the equation 

A(^,) = (4.15) 



is satisfied, with the additional condition that the first non-zero derivative of A{6) at Or is 
odd. 
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5 0+(2A^) symmetry 

We now address the question of finding a class of symmetries of the Hamiltonian (jH.lll 
which leads to an interpretation of the spin-spin correlation functions and the formula for 
the entropy of the entanglement ()3.1H) as averages over 0"'"(2A^), the group of orthogonal 
matrices of dimension 2N x 2N and determinant 1. We have seen that the expression of 
the von Neumann entropy in terms of an average over U(A^) is a direct a consequence of 
the invariance under translations of the Hamiltonian (jH.lll and of its geometrical isotropy. 
We now proceed in the same way as with U(A^) and try to infer how the structure of the 
kernel of the Haar measure of 0^(2A^) is reflected into the invariance properties of if^- 

Eigenvalues of orthogonal and symplectic matrices come in complex conjugate pairs, 
therefore 0+(2A^) has only N independent eigenvalues. When dealing with the classical 
compact groups, we shall adopt the convention of denoting by N the total number of eigen- 
values and by the number of independent ones. In general we shall denote an arbitrary 
group by 0(^1,^2) (^)- Each of the classical compact groups is identified by specific values of 
<^2)- This correspondence is described in appendix C; for 0~^{2N), (cji, £12) = (1/2, 1/2). 
Let F{U) be a class function on G(ai,a2){N), i-e- a symmetric function depending only on 
the eigenvalues of U. Furthermore, suppose that 



N 



F{u) = l[f{e,), 



(5.1) 



where f{6) is even and 27r-periodic. The averages discussed in appendix C can all be 
written as 



F{U) 




(5.2) 



j,k=l,...,N 



d^l ■ ■ ■ d^^. 



The quantity Q%^''^^\0i^92) is called the kernel oi the Haar measure and 



N 



det [/(^,)/(4)Q;:,.,) {0,M]^^,^,_^ = [\{f{0,)f{-e,) <,.,)(^i, . . . ,^;v), (5.3) 



where P^^^ o-2)(^i' • • • '^Jv) is the Haar measure ()C.25|) . The integral ()5.2p can always be 
expressed in terms of the independent eigenvalues; indeed in appendix C it is shown that 
it is always proportional to the integral 



N 



TT /"TT 



JO 




„...,^A.) d^i---d^ 



N 



(5.4) 



det(4^''"'^)j,fc=o,...,iv-i, 
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with a constant of proportionality that depends on the group and on the function f{6). In 
equation ()5.4|1 we have set g{6) = f{6)f{—6). Exphcit expressions for the matrix elements 
dp}''^^^ and the relations between the averages ()5.2|1 and ()5.4|1 for the various compact 
groups are reported in table 1, appendix C. In the rest of the paper we shall concern 
ourselves only with integrals of the form ()5.4|) . 

Let us now go back to 0^(2A^). In appendix C we show that 

"00 = 9o (5.5a) 
aoj = ajo = V2gj, j > 0, (5.5b) 
ajk = Qj-k + gj+k, j, k>0, (5.5c) 

where for simplicity we have dropped the superscript (o"i,(T2) and 

27r 



91 = 1^1 9{0)e-'''de (5.6) 



is the Ith Fourier coefficient of g{9). Matrices of the form {hj+k}j,k=o,...,N-i are called 
Hankel matrices, therefore the matrix {(yjk}j,k=o....,N-i is always the sum of a Toeplitz and 
a Hankel matrix. 

How can we infer from ()5.Hj) the structure of the matrices A and B that appear in the 
Hamiltonian ()3.ip ? After all, the geometry of is that of a discrete lattice while the 
kernel of 0''"(2A^) lives on the circle, its explicit form being 

Qo+(2iv)(0, ^) = S2N-i{<P -^) + S2N^i{<P + ^), <P,^e [0, vr), (5.7) 

where Sn{z) is the kernel ()4.2j) . In appendix C it is shown that the matrix elements ajk can 
be expressed as integral transforms involving a particular class of orthogonal polynomials, 
known as Jacobi polynomials. Furthermore, the kernel of Haar measure can be expressed 
in the form (see, e.g., [Hj, p. 24) 

Af-l 

Qo+{2N)i(P^^) = ^Pj{cos(j))pj{costlj), (5.8) 

j=0 

where Pj{x) is the jth Chebyshev polynomial of the first kind: 

1 [2 
Po(x) = —= and Pj(x) = \ — cos (j cos ^ x) , j > 0. (5.9) 
^/7^ V vr ^ ^ 



Formula ()5.8|1 leads to the following expression for the matrix elements that appear in the 
determinant of the left-hand side of ()5.3j) : 

N-l 

/(0)/(^)Qo+(2iv)(</>,^) = /(0)/(^) $^P,(cos</))p,(cos^). (5.10) 

j=0 

Then, if we compare the integral transforms ()C.35|) with the sum ()5.10|) . we note that the 
expressions are the same, but the role of the continuous and discrete variables is exchanged. 
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In other words, the functional form of the kernel ()5.7|) is complementary to the intrinsic 
structure of the matrix {ajk}j,k=o,...,N-i as a Toeplitz plus Hankel matrix. It is therefore 
natural to assume that the matrices A and B defining the quadratic form (jH.lj) should be 
the sum of Toeplitz plus Hankel matrices. As for U(A^), the periodic boundary conditions 
will impose on them a further structure which will turn out to be essential to our study. 
It is worth noting that the analysis presented in section ^ in terms of invariance under 
translations of the Haar measure of U{N) and of Ha is equivalent to the one discussed 
here; in the case of U{N) the orthogonal polynomials Pj (cos0) are replaced by the complex 
exponentials e'-''^/ ^/2n. 

The above considerations lead one to consider matrices A of the form 

Aji = a{j-l) + a{j + l), J, / = 0,...,M-1. (5.11) 

Because of the periodic boundary conditions a must be a function on Z/MZ, which must 
also be even in order for A to be symmetric. Clearly, a Hankel matrix cannot be antisym- 
metric, therefore 7 must be zero: the Hamiltonian (jH.lj) must be isotropic. A brief look to 
table 1 in appendix C shows that this is a necessary condition for all the other compact 
groups too. Since A is a real symmetric matrix, its eigenvalues A^ are real and therefore 

V'fe = ^0fc = signAfc^fc, (5.12) 

where the 4>j^s are the eigenvectors of A. 

We now need to diagonalize A; as for the unitary group, this can be done explicitly. 
Because 

Aji = AjM-h (5.13) 
any odd function on Z/MZ will be in the kernel of A. Therefore, 

27r/ 




,, = y — sm^-j, k = —, Z = l,...,[(M-l)/2], (5.14) 

are a set of independent eigenvectors with eigenvalue zero whose multiplicity is at least 
[(M — l)/2], where [■] denotes the integer part. 

The eigenvectors with non-zero eigenvalues can be found by exploiting the symmetries 
of the matrix ()5.1H1 . For simplicity, we assume that the non-zero eigenvalues are non- 
degenerate. An immediate consequence of ()5.1H|1 and of the condition A^ 7^ is that any 
eigenvector Ck{j) must be an even function on Z/MZ; thus, we can always write 

Cfc(j) = efc(j) + efc(-j), J e Z/MZ. (5.15) 

Furthermore, using the periodicity of a(j), it is easy to show that if c(j) is an eigenvector 
of A, than Ck{j + p) + Ck{j — p), where p is an arbitrary integer, is also an eigenvector 
corresponding to the same eigenvalue A^. It follows that 



CkU +p) + Ck{j -p)gc e',,{p)ck{j). 



(5.16) 
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Since Cfc(j) is even, the role of j and p can be interchanged, and since they are both 
arbitrary, we can choose e'^(j) = efc(j). Then it follows from equations ()5.15j) and ()5.1(i|l 
that 

efc(j +p) = efc(j)efc(p), (5.17) 

for an appropriate choice of the constant of proportionality in ()5.16p . Thus, the efc(j) are 
additive-multiplicative functions. Because they are periodic too, it must be that 

ek{M)=ek{l)'' = 1. (5.18) 

Therefore the efc(j)s are roots of unity: 

p^'^^ 27r/ 

efc(j) = ^, k=—, J,/ = 0,...M- 1. (5.19) 



As immediate consequence of ()5.15|) the remaining normalized eigenvectors of matrix ()5.11|) 

= A' '^^i = \li < A; < TT for M odd 

= 7m ' = yn^^^ < A; < vr, (f)„j = for M even, 

where k = 27r//M. The corresponding eigenvalues can be obtained by direct substitution: 

^ ^ f2a(0)+4Ei=r'^/'«(j)cosA;j if M is odd 

^ [2 [a{0) + {-ly a{M / 2)]+ 4 J2fJi'^aU) cos k J if M is even. 

In appendix B we show that there exists a canonical transformation of the Fermi oper- 
ators bj that diagonalizes (jSH}. Using the same notation as in appendix B, let us denote 
by rjk the Fermi operators in term of which Ha is diagonal. The fact that approximately 
half of the A^s are zero means that the corresponding r^^s do not appear in Ha- In other 
words. Ha is isomorphic to a system with half the number of degrees of freedom. This is 
not surprising; it is a reflection of the fact that only half of the eigenvalues of a matrix in 
0+(2A^) are independent. In the same way as statistical properties of orthogonal and sym- 
plectic matrices are computed only in terms of the independent eigenvalues, so the extra 
degrees of freedom in Ha can be ignored. The matrix ()H.8|) can therefore be determined 
from the eigenvectors ()5.20|) . 

Following the same steps as for U(A^) we fix our attention on the subsystem P composed 
of the first consecutive oscillators and let M ^ oo. The eigenvalues 1)5.2111 converge to 
the even function 

oo 

A{e) = Ao + 2 ^ Aj cos ej, (5.22) 

i=i 

where Aj = 2a{j). Finally, by substituting the vectors ()5.20|) into ()3.8|) and taking the 
limit M — i> cxD, we obtain 

2 Afc 

(7iv)i7 = — > 7": — ; coskj cosklAk > (5.23) 

71 ^ \Ak\ M->oo 
k=0 ' ' 



Random Matrix Theory and Entanglement in Quantum Spin Chains 



19 



where the ajk are precisely those of equation ()5.5|) with symbol 



9(0) 



(5.24) 



\m\' 



We can then define an XX model with orthogonal symmetry by choosing 



a{j) = aMj) + 5,{~j)]-25oij). 



(5.25) 



It is important to notice that in order for T^r [g] to be the sum of a Toeplitz and a Hankel 
matrix the subchain P must be made of the first N sequential oscillators: we cannot shift 
the subsystem P because Ha is not translation invariant. This property was to be expected 
because the kernel ()5.7|) is not invariant under translations. In other words, the origin of 
the lattice defining the spin chain is a privileged point, in the same way as the point 1 on 
the unit circle is a symmetry point for the spectra of orthogonal and symplectic matrices. 
A consequence of the absence of translational invariance is that the determinant expressing 
the spin-spin correlations ()H.7j) is an average over 0"'"(2A^) only if one of the spins is the 
first in the chain. 

It turns out that the structure of the matrix T^lg] so obtained has important conse- 
quences for the formula ()3.1H) . First, since A(— ^) = A(^), Tj\f[g] is symmetric, therefore 
the correlation matrix factorizes as in ()2.23|) . Thus, the characteristic polynomial 



in the integral (jH.llj) is an average over 0^(2A^). We compute this integral in section|Hl The 
symbol A{6)/\A{6) \ is the same as the one discussed in section |3] for the case when 7 = 0: 
it is a piecewise continuous function that takes the values 1 and —1 and whose jumps are 
located at the points 6r which are solutions of equation ()4.15|) . In section 0] we have seen 
that the eigenvalues of the corresponding Toeplitz matrix are in the interval [—1, 1]; similar 
arguments lead to the same conclusion for the eigenvalues of a matrix which is the sum of 
a Toeplitz and a Hankel matrix with the same symbol. Therefore, formula p.lip gives the 
entropy of the subchain P and can be expressed in term of an average over 0^(2A^). 



The treatment of these two groups turns out to be the same - see (IC.49d|) . The arguments 
are analogous to those presented for 0"*'(2A^). The elements in Sp(2A^) are 2Nx2N unitary 
matrices U such that 



where I is the N x N identity matrix. The number of independent eigenvalues in both 
0~(2A^ -|- 2) and Sp(2A^) is N. Without loss of generality, we shall concentrate only on 
Sp(2A^). The kernel of the Haar measure 



DN[g]{X) = det{XI-T!,[g]) 



(5.26) 



6 Sp(2iV) and 0-{2N + 2) symmetry 




(6.1) 



<5sp(2iV)(0> i') = S2N+li(p - ^) - S2N+l{(j) + ^), (j),^ e [0, Tt) 



(6.2) 
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and the matrix elements appearing in the determinant ()5.4|) are 

Oijk = Qj-k — gj+k+2- (6.3) 

As for 0^(2A^) the structure of formulae ()6.2|1 and ()6.3|1 is incompatible with 7 7^ 0. 
The choice of the matrix A is 

'Ajk = a{j -k)-a{j + k + 2), (6.4) 

where a is an even function on 'L/M'L. The diagonalization of the matrix ()fi.4|l is analogous 
to the one of ()5.1H) . Therefore, we just present the results. The eigenvectors that span 
the kernel of A are 

[Y 27r/ 
0fei = y^cosfc(j + l), k=—, / = 0,...,[M/2], (6.5a) 

while those corresponding to 7^ are 

T _ . . 27r/ 




= y-sin^j + l), fc = — , / = l,...,[(M-l)/2]. (6.5b) 

As in the case of 0"'"(2A^) the Hamiltonian ()3.H1 is isomorphic to a system with half 
the number of degrees of freedom. The eigenvectors associated to the relevant degrees 
of freedom are those in ()6.5bjl . Similarly, the eigenvalues can be computed by direct 
substitution; they turn out to be given by formula ()5.21|) . By fixing the number of 
oscillators in the sub chain P and letting M — 00, the matrix T/v converges to 

(Tm),i = - ^ sin k{j + 1) sin k{l + 1) Afc 

u—n I 



fc=0 



(6.6) 



27r 7o |A(^)| 

which are the integral transforms ()(I39jl : A(^) is the same real and even function as 
in ()5.22|) . therefore Tfq[g\ is symmetric. An immediate consequence is that the entropy 
formula ()3.11|) can be expressed in terms of an average over Sp(2A^). 



7 0^(2iV + 1) symmetry 

The treatments of 0''"(2A^ + 1) and of 0^(2A^ + 1) follow a similar pattern, indeed averages 
over these two groups are intertwined by equations (jC.49bjl and ()(I49cj) . The kernels of 
the Haar measures are 

Qo±(2^+i)(0, V') = ^2^(0 ^2jv(0 + V^), 0, e [0, tt), (7.1) 

and the matrix elements in the average ()5.4|) are 

otjk = gj-k T dj+k+i, (7.2) 
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where in equations ()7.1|) and ()7.2j) the minus sign refers to 0"^(2A^ + 1) and the plus sign 
to 0~(2iV + 1). Consequently, the choices of the matrix A compatible with ()7.H) and ()7.2|1 
are 

Ajk = a{j -k)Taij + k + l), (7.3) 

where a is an even function on Z/MZ. As with the groups treated previously, equa- 
tions ()7.1|1 and ()7.2|1 are incompatible with 7 7^ 0. 

The matrix ()7.3|1 can be diagonalized adopting the same techniques used for the other 
groups. 




M 
2 



sm 




M 



cos 



2j + l 
2 

2j + l 



k ■ 
k 



27ll 

2txI 



/ = l,...,[M/2], 

/ = 0,...,[(M-l)/2]. 



(7.4a) 
(7.4b) 



These are the eigenvectors of the matrix ()7.3|) for both choice of signs; however, the func- 
tions ()7.4b|) are in the kernel of A when the sign between the two terms in ()7.3|) is minus, 
i.e. for 0"'"(2A^+ 1) symmetry, while their eigenvalues are not zero when the sign is plus, i.e. 
for 0~(2A^ + 1) symmetry. For the eigenvectors ()7.4a|l the role is reversed. The eigenvalues 
are given by formula 1)5.2111 for these groups too. It is now straightforward to determine 
the matrices T/v for both groups: 



{Tn) 



2 \ ^ Afc . 

71 ^ Afc 

fc=0 ' 



k 



2j + l 



sm 



k 



21 + 1 



1 



- m 

M^oo 271 Jo \A{e)\ 



(7.5a) 



(e"iO-Oe _ e-Ki+'+DS) 



for 0+(2iV+ 1) and 



2^ A, 

- / "TTT COS 



2j + l 



cos 



21 + 1 



Ak 



M^oo 27r 7o |A(^)I 



(7.5b) 



for 0-(2A^ + 1). The function A(6') is even and the symbol g{d) = A{e)/\A{9)\ is the 
same function analyzed in connection with the other groups. The matrix T]^[g] is real and 
symmetric and therefore the formula ()3.11|1 for the entropy of the entanglement has an 
interpretation as average over 0"'"(2A^ + 1) or over 0~(2A^ + 1). 



8 Generalizations of the Fisher-Hartwig formula and 
the computation of entanglement 

The computation of Toeplitz determinants, and in particular of their asymptotics, is im- 
portant in many branches of Physics. The first and most famous application goes back to 
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1946 and is due to Osanger, who showed that the diagonal spin-spin correlations in the 
classical two-dimensional Ising model can be expressed in terms of Toeplitz determinants. 
It turns out that the behaviour of the leading order term as the dimension of the matrix 
tends to infinity changes radically when the symbol has discontinuities or zeros. Indeed, 
phase transitions in quantum and classical lattice systems often appear as changes in the 
analytic properties of symbols in Toeplitz determinants. 

When the symbol g{6) is a continuous function on the unit circle and the c^s are the 
Fourier coefficients of log g{9), Szego's theorem jEj states that 



lndetrAr[5(] = coA^ + ^fccfcC_fc + o(l), ^ oo, (8.1) 

k=l 

provided that the series YlT=-oo\^k\ and X^fcl-ool^ll'^fcl^ converge. If g{6) has zeros or 
discontinuities, then it can always be reduced to the form 



g{e) = mH^aAS - dr)tpXd - dr), (8.2) 



r=l 



where (f) is smooth, has winding number zero and 

tp{e) = exp [-ip (n-e)], 0<9<2n, /3 ^ Z (8.3a) 
ti^(0) = (2-2cos^)", Rea>-^. (8.3b) 

Note that L represents the number of zeros/discontinuities in the interval [0, 27r). Fisher 
and Hartwig ffl conjectured that 



lnDM[g] = coN+ \^^al-(3^jlnN + lnE + o{l), N ^ oo, (8.4) 

where now the c^s are the Fourier coefficients of In 4>{6). Basor determined the constant 
E: 



E = exp 5^ kcc, n (^''0 "'"^^ <P- i^-''') ~' ' 

\k=l J r=l 

X n (1 - exp [i {e, - g,)])-K+/^-)("--/^-) -Q G (1 + a. + ^) G (1 + ar - (3r) 



G (1 + 2ar) 

l<r=^s<L r=l ^ ^' 



1.5) 
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where G{z) is the Barnes G-function 



? and 



oo 



oo 



ln0+(t) = Y^c.t^, ln0_(t) = Y,c^,t-K 



(8.6) 




The Fisher-Hartwig conjecture has been proved for |Rear-| < 1/2 and |Re/3r| < 1/2 jUj 
and for other specific values of a^, /^r and L. (The reader is referred to [TH], pp. 469- 
474, for a complete discussion.) The most important difference between equations ()8.1|) 
and ()8.4|) is the extra term proportional to In N in the Fisher-Hartwig formula. 

Formula ()8.4|) was used by Jin and Korepin |1] to compute the entropy of the entangle- 
ment for the XX model, that is for the Hamiltonian 1)2.11) with 7 = 0. Their computation 
can be easily generalized to all Hamiltonians of the form p.ip that are invariant under 
translations and isotropic, for which, therefore, the formula ()3.11|1 becomes an average 
over U(iV). However, before entering in the details of the computation, we discuss the 
generalization of formula ()8.4j) to determinants of combinations of Toeplitz and Hankel 
matrices, or, more specifically, determinants that can be interpreted as averages over the 
other classical compact groups. 

Basor and Ehrhardt proved a generalization of the Fisher-Hartwig formula to deter- 
minants of matrices of the type 



i.e. determinants that are averages over 0-(2A^ + l). Using equations (l(J.49bD and (l(J.49c;i . 
their formulae can be applied to averages over 0^{2N -|- 1) by a change of variable in 
the integral ()5.3|1 . In a recent paper Forrester and Frankel p] extended the results of 
Basor and Ehrhardt conjecturally to averages over Sp(2A^), 0^(2A^) and 0~(2A^ + 2). 
These formulae provide a very efficient means to compute the integral ()3.11)1 for chains of 
fermionic oscillators whose symmetries are associated to one of these groups. 

From the discussion of the previous sections it has emerged that we are only interested 
in the symbol g{9) = A(^)/|A(^^)| when A(^) is real. This symbol takes only two values: 1 
and —1. Its discontinuities are located at the points where equation ()4.15p has solutions. 
If such an equation has no solutions, then g{6) is a constant and Tj^[g] = ±/ for all the 
compact groups. It follows from equation ()2.2U|) that E-p = 0. In physical terms this means 
that the Hamiltonian Ha is away from the critical point or, equivalently, the magnetic field 
is so strong that all the spins are aligned, thus there are no correlations and correspondingly 
the entanglement must be zero. 

The Fisher-Hartwig formula extended to all the compact groups is fairly complicated 
when expressed for a general symbol of the form ()8.2)1 . For simplicity, we shall report 

^The definition of the Barnes G-function is 



9j~k + 9j+k+U j, = 0, . . . , iV - 1, 



(8.7) 



G(z) = (2^)^/2 e-[.(.+i)+7E.^]/2 -Q ^(^^ ^ 



-z+z^/{2n) 



n=l 



where 7e is Euler's constant. It generahzes the Gamma function, in the sense that it obeys G(z + 1) = 

r(z)G(z). 
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only the expression for the case that concerns us directly, i.e. g{6) is even and has only 
discontinuities; for the general case we refer to [F. Equation ()8.2|1 becomes 

R 

g{e) = mUt^rio - er)t.pXo + Or). (s.s) 

r=l 

where now all the discontinuities 9^ lie in the interval [0, vr) and therefore L = 2R] we 
exclude the case 6*^ = 0, vr and still require |Re/5r| < 1/2. We have 

lnD;vM(A)K,.,) = Nco - (jZ A') InN 

^'=' ^ (8.9) 

^ oo oo ^ I 

+ 2 5Z + ^2fc-i + In ^K,<x2) + In E + o(l), AT ^ oo, 

k=\ k=l 

where 

R 

<T2) = 1^(1 + e'^'-)'^''l"^/^)^'-(l + Q-''Qry{ai~l/2)fir(^l ^ gi9.^(a2 + l/2)/3,j^-|^ _^ ^_i9,^-(a2 + l/2)/3. 



r=l 



exp <^ - ^ Cfc [(71 - 1/2 + (-1)^= (a2 + 1/2)] I (8.10a) 



fc=i 



JL , nil- e-i«r■^/3./2/'1 _ ie,N-/3,/2 



r=l 

2ft./3a R 



r=l 



l<r<s<_R 



Note that the main differences between the various groups appear only in the 0(1) 
terms of ()8.9|) : when an even symbol with discontinuities is averaged over a compact group 
the term linear in in formulae ()8.4j) and ()8.9|1 is the same in every case, and the term 
logarithmic in A^ has just an extra factor of two in front of it for U(A^). This factor is due 
to the fact that averages over U(A^) are computed by integrating over [0, 27r)^, while for 
the other compact groups the domain of integration is [0, tt)^, therefore, in the latter case 
the singularities located at —9^ do not contribute. 

The representation (j8.8j) of g{9) = A{9)/\A{9) \ is given by the following choices of - 
and Pr- 



0(A) = (A + l)(^^j , (8.11a) 

A.(A) = i-iym, m = ^ In , (8.11b) 

with —IT < arg[(A + 1)/(A — 1)] < vr. Therefore, we have |Re/3r| < 1/2 on the contour of 
integration c(e, 5) and we can apply formula ()8.9|) . The leading order asymptotic of the 
entropy (jH.llj] is then given by 

Ep = hN -2'"^Rl2\nN + 0(1), N ^ oo, (8.12) 
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where R is the number of discontinuities in the interval [0, tt) and 

fl if the average is over U(A^) 
^ ^ ' (8.13) 

if the average is over the other compact groups. 



Ji and I2 are the integrals 



/i = lim lim — e(l + e, A)^:^ dA = 0, (8.14a) 
/2=limlim — / e(l + e,A)/5(A)/?'(A)dA = -l/(61n2). (8.14b) 

,^0+ 5-^0+ TTl ,5) 

The first integral can be computed straightforwardly using the residue theorem; the second 
was computed in j4jj (in this case the integrand is a multivalued function inside the contour) . 
Finally, we have 

210G R 

Ep log2 A^, ^ 00. (8.15) 

6 

The asymptotic relation ()8.15|1 represents one of our main results. As mentioned in 
section |21 the logarithmic growth of the entanglement is a general consequence of the 
fact that in one dimensional systems near quantum phase transitions the entropy is a 
logarithmic function of the size of the system [SJ IH] . These thermodynamic arguments give 
the value of the coefficient in front of the logarithm to be one-third of the central charge of 
the associated Virasoro algebra. Our expression therefore leads to an explicit formula for 
the central charge, which depends in a non-trivial way on the geometry of the Hamiltonian. 
The factor 2"'^: is universal, depending only on the symmetries determining the classical 
compact group to be averaged over. The factor R corresponds to the number of solutions 

of Km . 

We now proceed to compute the next-to-leading-order term. This can be determined 
by integrating with respect to A the terms independent of A^ in the Fisher-Hartwig formu- 
lae (Q and (jHini)- We begin with U(A^). 

Since the symbol is even, using ()8.8p and ()8.11ap . the constant ()8.5p becomes 



i?(A) = J]|l-e'2^'-|-2'3W' n 



1 _ e'ier-9s) ^-D^^-^-'^m' 



X [G(1+/?(A))G(1-/3(A))]^ 



r=l l<r<s<B. 

i2R 



(8.16) 



where only the independent discontinuities, located in the interval [0,7r), are taken into 
account. The logarithmic derivative of ()8.1(i|l is the sum of two terms. The first one is 

4/?(A)/5'(A) ("2 Y: (-l)^^^^^l- -ElHl-e'^^l) (8.17) 



l<r<s<_R 

The second term is more delicate; we have 



r=l 



G (1 + /3(A)) G (1 - m) = e-^W'(i+^-) f] (l - ^) ^ e 

n=l ^ ^ 
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where 7e is Euler's constant. The logarithmic derivative of the right-hand side of ()8.18|) is 

-2/?(A)/?'(A)[l + 7E + T(A)], (8.19) 



where 



Let us define 



n=l ^ ' 



1 / ^ 

ir = l + 7E + ^ ^ln|l-e^2^>.|_2 ^ (-!)('■+«) In 



R 

and 



,r=l l<r<s<_R 



1 _ f^ASr-es) 



1 _ eKfr+^s) 



(8.21) 



/3 = hm hm — / e(l + e,A)/3(A)/3'(A)T(A)dA. (8.22) 

e^0+ 5^0+ TTl ,5) 

This integral was evaluated in |^ where it was found that /a = 0.0221603.... 

Combining equations ()8.16|) . ()8.17|) and ()8.19|) . we obtain that the 0(1) contribution to 
the entropy of entanglement is 

Cu(Ar) = 2RhK - 2Rh = ^ {K - 6/3 In 2) . (8.23) 

When R = 1 this equation reduces to the result of Jin and Korepin for the XX model. 

Let us now consider the other compact groups. From equations ()8.1()|1 it is clear that the 
next-to-leading-order term is composed of two parts, one, common to all groups, coming 
from equation ()8.10b|) and the other, depending on the particular choice of the group, 
given by ()8.10a|) . By taking the logarithmic derivative of equation ()8.10b|) one immediately 
realizes that only two integrals contribute to this term: the first one is I2', the second one 
is 

J4 = -L/ e(l + e,A)/3'(A)dA. (8.24) 

'c(e,(5) 



27ri 



This integral can be evaluated using the residue theorem to obtain I4 = 0. Therefore, by 
proceeding as for U(A^) we can determine the contribution to the sub-leading term coming 
from equation ()8.10b|) : 

Cg = ^(i^-6/3ln2 + ln2). (8.25) 

This expression differs from that for Cu(Ar) by a factor 1/2, the origin of which is the same 
as in the leading order term, and an additional term equal to R/6. 

We are left to determine the contribution coming from equation ()8.10ap . It is evident 
that, up to a constant depending only on the ^^s, by taking the logarithmic derivative of 
the right-hand-side we are left just with integrals of the type fj8.24|) . Hence, the next-to- 
leading-order term for the groups Sp(2A^) and O^(A^) is simply ()8.25|1 . 
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9 Painleve VI and gap probability generating func- 
tions 

The averages we have discussed turn out to be related to solutions of integrable second order 
ODEs of Painleve type. There exist six Painleve equations; any second order differential 
equation free of moveable essential singularities of the form 

y" = R{y',y,t), (9.1) 

where i? is a rational function, can be reduced to one of them. 

Let us consider a generating function g{6) = A{6)/\A{6)\ with only two discontinuities 
at 6i and —6i and set (p = 26i. We then look at the characteristic polynomial ()3.12j) when 
g{6) is averaged over U{N), which can be written as the integral 

[ZTV) l\. Jo Jo .^^ L J ^^.^^^^ 

= (A + 1)^ E^^^{Ii^_^^^y, ^), 

(9.2) 

where we have set ^ = 2/(A + l). The function i?^^^(/[^_0^^); C,) is known in random matrix 
theory as the generating function of the gap probabilities. Here CUE stands for Circular 
Unitary Ensemble, which, as already noted, denotes the probability space given by U(A^) 
equipped with the Haar measure. Differentiating one obtains 



En = [-1) 



dC 



(9.3) 

€=1 



where E^^^{n; /[^„0 ,r)) is the probability that the interval /[^r-^^Tr) contains n eigenvalues. 
Making the substitution 

e'^ = , z = tan - (9.4) 

1 — ia; 2 



equation ()9.2|) becomes 



N 



^^-f[s,oo)(-^i) 



X Y\ \xj — Xkf dxi ■ ■ ■ dxN, 

^<j<k<N 

where s = cot (0/2). The integral ()9.5|) is the gap probability generating function of the 
Cauchy ensemble (CylJE) for the interval [s, oo). 

If we let the exponent in the denominator in the integrand ()9.5j) vary, we can define 



c\N^ poo poo ^ 



'j=i i^+^'iy (9.6) 

X Y\ \xj — Xkf dxi ■ ■ ■ dxN. 
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The connection between the gap probabihty generating function of U(A^) and the theory 
of Painleve equations is mediated by the function 

a{s) = (1 + s^)^ In a + N,^, (9.7) 
which is a solution of the equation 

(1 + s^fiay + 4(1 + s^)iaf - %sa{o'f + 4a2(a' - a^) 

+8a2saa' + 4[A^(A^ + 2a) - a^s^\ [a'f = 0. (9.8) 

The above equation is known in the hterature as the cr-form of the Painleve VI equation. 

Solutions of equation ()9.8|1 . and therefore the determinant ()9.2j) . obey a recurrence 
relation which allows one to determine an exact formula for them for each value of A^. 
Therefore, at least in principle, it is possible to compute the entropy of the entangle- 
ment ()3.11|) exactly for each A^. In terms of of the gap probability generating function, 
these recurrence relations assume the following form |19j : 

rnCUE pCUE 

^ ^ _ 2 g^ 

with initial conditions 

271 

In turn xn obeys the recurrence relation 

1 — 



2xnXn-1 - 2 cos ^ = — [{N + l)xAr+i - (A^ - l)xAr_i] 

2 Xn 



[NXN -{N- 2)XN^2] 



(9.11) 



with initial conditions 

^ sin ~ 

x_i = 0, Xq = 1 and Xi = 1 — . (9.12) 

These recurrence relations can be used to compute higher order terms in the asymptotics 
of E^^^ as ^ oo, and therefore of the entropy of the entanglement. Substituting the 
Fisher- Hart wig formula ()8.4j) into ()9.9|1 gives 

A,_oo. (9.13) 

where /?(A) was defined in ()8.11aj) . This suggests that xn has an asymptotic expansion of 
the form 
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Inserting this expansion into ()9.11|) it is possible to compute recursively the coefficients 
Cj{X); for the second order term we obtain 

C2(A)=2^/3|/3(A)|l (9.15) 

This coefficient determines the contribution to the integral (jH.llll given by equation ()8.2;^j) . 
Higher order terms can be computed in a similar way. 

When the average is over the other compact groups, the determinant ()3.12|1 can still be 
interpreted as gap probability generating function for the respective group and its logarithm 
is still a solution of a differential equation related to the Painleve VI equation. However, 
although recurrence formulae analogous to ()9.9|) exist, at each step the values of (cri,(T2) 
that label the integral ()5.4|) change and in general do not even identify one of the classical 
compact groups (see, e.g., [20]). 

10 Conclusions 

We have investigated the entanglement of formation of the ground state of the general 
class of quantum spin chains related to quadratic Hamiltonians of fermionic oscillators 
partitioned into two contiguous subchains. The number of oscillators in the ffist subsystem 
is and we let the total size of the system grow to infinity. We have discovered that 
for certain Hamiltonians the measure of entanglement, which in these circumstances is 
given by the von Neumann entropy of the first subchain, can be expressed in term of 
an average over one of the classical compact groups U(A^), Sp(2A^) and O^(A^). Indeed, 
we show that there exists a one-to-one correspondence between the symmetries of the 
fermionic chain and the functional form of the Haar measures of the classical compact 
groups. Using generalizations of the Fisher-Hartwig conjecture it is possible to compute 
asymptotic formulae for such averages in the limit N —>■ oo. The entanglement is either 
zero, away from critical point, or grows logarithmically with in the proximity of quantum 
phase transitions. Generalizations of the Fisher-Hartwig formula allow one to determine 
the constant in front of the logarithm explicitly, and the next-to-leading-order term in the 
asymptotics. Furthermore, these averages turn out to be related to solutions of Painleve 
equations. 

The fact that one can compute the leading order terms in the asymptotics of the en- 
tanglement of formation of the ground state of such a significant class of systems suggests 
that the connection between lattice models and random matrices may be deeper than be- 
ing simply a calculational device. For example, the diagonal spin-spin correlations of the 
two-dimensional classical Ising model are Toeplitz determinants with symbols analogous to 
the one associated to the XY model; it is likely that a similar association between classical 
compact groups and symmetries of the Hamiltonian exists also for classical lattice mod- 
els. If it does, random-matrix techniques could then be used to deduce thermodynamic 
quantities like critical exponents. After all, two-dimensional classical spin chains are math- 
ematically equivalent to one- dimensional quantum spin chains. Furthermore, these type 
of random matrix averages already appear in the calculation of the ground state density 
matrices for an impenetrable Bose gas in an interval of finite length [Oj. 
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It is likely that these applications of group averages will prompt further studies of 
the analytical properties of the determinants and spectra of combination of Toeplitz and 
Hankel matrices, the investigation of which has started only recently. 
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Appendix A. The density matrix of a subchain 



Let {l^/'j)} be a basis of the Hilbert space 7i of a system composed of two parts, P and Q, 
so that 7i = Tip (8) Tiq. The density matrix of a statistical ensemble expressed in the basis 
{IV'j)} is a positive Hermitian matrix given by 



jk 



(A.l) 



with the condition Trpq ppq = 1. Let us introduce the operators S{j, k) and S'(j, k) defined 
by the relations 



s{i, k) = \ipj){iJk\ 

S{i,k)S{k,l) = 5ji\il)j){il)j\ and S{i,k)S{k,l) = 5ji\ipj){il)j\ 
(In this formula repeated indices are not summed over.) Clearly, we have 

Cjk = TrpQ [ppQS{k,j)] . 



(A.2a) 
(A.2b) 

(A.3) 



Let us now suppose that the Hamiltonian of our physical system is ()3.1|) and that 
the subsystem P is composed of the first oscillators. Then a set of operators S{j, k) 
for the subchain P can be generated by products of the type H^Li^'j, where Gj can be 
any of the operators {cj,CpC^jCj,CjC^j} and the Cjs are Fermi operators that span Hp] it is 
straightforward to check that S{k,j) = We then have 



Pp = '^^P 

All the S{1, k) 

= E Trp 

All the S(l, k) 

= '^^PQ 

All the S{1, k) 




(A.4) 
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Since ppq = |^'g)(^'g|, this expression simply reduces to 

N N t 



All the S{1, k) 



N 



The correlation functions in the above sum can be computed using Wick's theorem ()3.7j) . 
Finally, if the correlations of the CjS are given by ()2.8|) . we immediately obtain formula ()2.18|) . 

Appendix B. The correlation matrix Cm 

The purpose of this appendix is to provide an explicit derivation of the expectation values 

(*g|mjmfc|*g) (B.6) 

when the dynamics is determined by the Hamiltonian ()3.1|) . 

First, we need to diagonalize ifa, which is achieved by finding a linear transformation 
of the operators hj of the form 



M-1 

r]k = 

j=0 



5^ [gkjb, + hkjb^^ , (B.7) 

j=0 

such that the Hamiltonian ()3.1|) becomes 

Af-l 

H^=Y,\^k\vivk + C, (B.8) 

k=0 

where the coefficients Qkj and h^j are real, the rj^s are Fermi operators and C is a constant. 
We use the notation \Ak\ because it is convenient for the computations carried out in 
section Hlto allow the coefficients of the number operators rjlrjk to be the absolute values of 
the complex numbers A^. The quadratic form ()3.H) can be transformed into ()B.8|1 by ()B.7|) 
if the system of equations 

[7]k,H^]-\Ak\Vk = 0, k = 0,...,M-l (B.9) 

has a solution. Substituting ()3.1|) and ()B.7|) into ()B.9|) we obtain the eigenvalue equations 

M-1 

\^k\gkj = ^ {gkiAj - hkiBij) , (B.lOa) 

(=0 

M-1 

\^k\hkj = E {gkiBij - hkiAij) , (B.lOb) 

1=0 

where A = a A — 21 and B = a'jB. These equations can be simplified by setting 

(pkj = gkj + hkj (B.lla) 
i^kj = gkj - hkj, (B.llb) 
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in terms of which the equations ()B.10|) become 

(Z + 5)0, = lAfclV', (B.12a) 

(A-B)Vfe= |Afc|0,. (B.12b) 

Combining these two expressions, we obtain 

(A-B)(A + B)^, = \A,\'<I>, (B.13a) 

(A + B)(A-B)^P, = lAfclVfc. (B.13b) 

When Ak ^ 0, 0^ and \Ak\ can be determined by solving the eigenvalue equation ()B.13ajl . 
then -j/?^ can be computed using ()B.12a|) . Alternatively, one can solve equation ()B.13b|) 
and then obtain (pj^ from ()B.12b|) . When A/^ = 0, and i/?jr. differ at most by a sign and 
can be deduced directly either from (jB.12|l or from (jB.13j) . 

Since A and B are real, the matrices [A — B){A + B) and {A + B){A — B) are symmetric 
and positive, which guarantees that all of their eigenvalues are positive. Furthermore, the 
0^s and i/'fcS can be chosen to be real and orthonormal. As a consequence the coefficients 
Qkj and /ifcj obey the constraints 

A/-1 

X] idkjgn + hkjhki) = 5ji, (B.14a) 

fc=0 
M~l 

(gkjhki + hkjQki) = 0, (B.14b) 

k=0 

which are necessary and sufficient conditions for the r^^s to be Fermi operators. 

The constant in equation ()B.8|) can be computed by taking the trace of Ha using the 
two expressions ()3.H) and ()B.8j) : 

A/-1 A/-1 

TiHa = 2^'-' J2 ("^'^^ - 2) = 2"'-' E l^'^l + 2'''^- (B-15) 

k=0 k=0 

Therefore, we have 

^ M-l 

C = -^Mfc,-2-|A,|). (B.16) 

k=0 

We are now in a position to compute the contraction pair (jEH). Substituting (IrTT|) 
into ()B.7|) we have 

M-l 

r/fc = - ^ (0fci"^2j+i - i^kjm2j) . (B.17) 

j=0 

Since the 4>f^s and V'fcS are two sets of real and orthogonal vectors, ()B.17|) can be inverted 
to give 

A/-1 

m2j = i i^kj (vk - vl) (B.18a) 

fc=0 
A/-1 

m2j+i = 4>kj (vk + vl) ■ (B.lSb) 

k=0 
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Since the vacuum state of the operators rjk coincides with |^'g), the expectation values ()B.6|) 
are easily computed from the expressions ()B.18jl . We have 



M-l 



(*g| m2jm2k |*g) = ^ i'lj^ik = 5jk, 



1=0 

M-l 



and 



(*g| m2j+im2k+i |*g) = ^ (l)ij(f>ik = Sjk 

1=0 



M-l 

(*g| m2jm2k+i |*g) = i ipij(l)ik, 

1=0 
M-l 



(*g|m2j+im2fc |*g) = 
Finally, by introducing the real M x M matrix 



-i ^ i^ik4>ij- 

1=0 



M-l 



M)jk = ^l^k^ 
1=0 



j,k = 0,...,M-l 



and combining the expressions ()B.19|) we obtain 

(*g| mjmk |*g) = Sjk + i{CM)jk, 
where the matrix Cm has the block structure 



C 



( C\\ C\2 

C21 C22 



M 



CiM \ 

C2M 



\Cmi Cm2 ■■■ CmmJ 



with 



C, 



jk 



(B.19a) 



(B.19b) 



(B.20) 



(B.21^ 



(B.22a) 



(B.22b) 



{TM)jk 

-{TM)kj 

We call Cm the correlation matrix. It is worth noting that because of the definition ()B.20|) . 
the matrix Tm contains all of the physical information relating to the ground state of Ha. 



Appendix C. Averages over the classical compact groups 

Let f{6) be a 27r-periodic even function and let G(^„^^^^){N) be one of the classical compact 
groups 0+(2iV), 0-(2A^ + 2), Sp(2iV), 0+(2iV + l) and 0-(2A^+l), where the superscript 
± denotes the connected component of the orthogonal group with positive and negative 
determinant respectively. The integer N denotes the total number of eigenvalues, while N 
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denotes the number of independent ones. For any U G G(o-^^o-2)(^)) let ^{U) be the class 
function defined by 

TV 

F{U) = l[f{e,). (C.23) 

i=i 

In this appendix we want to compute averages of the type 

(FiU)) =(n/(%)) ■ (C.24) 

The Haar measure of G(^ai,a2){^) expressed in terms of the independent eigenvalues is 
given by (see, e.g., |2I], pp. 218 and 224) 

1 ^ 

'=1 (C.25) 

X JJ^ (cOS^j - COS^fc)^ 

l<i<'c<A^ 



where Z^^''^^^ is a normalization constant whose specific value will not be relevant to what 
follows. The parameters (ai, cxa) = (-1/2, -1/2), (1/2, 1/2), (-1/2, 1/2), (1/2, -1/2) refer 
to 0+(2A^), Sp(2A^), 0+(2A^+l) and 0~(2iV + l) respectively; we will consider 0-(2A^ + 2) 
separately. 

Averages of the type ()C.24|) can be evaluated by making the substitutions 

Xj = cos9j, (C.26) 
which reduces ()(I24jl to the computation of the integral 



1 rW ^ \ ^ 

^^(^i'^2) = -(^ / ■■■/ n^(^°^"'^) 11(1 + 

-^-1 "^-1 \j=l / i=l 

X Y\ {xj — Xk)"^ dxi ■ ■ ■ dxN, 

^<j<k<N 



(C.27) 



where g{6) = f{9)f{—9). For ai,a2 > —1 the integral ()(127jl can be evaluated using or- 
thogonal polynomial techniques; the details of the computation can be found, for example, 
in sections 2.1 and 2.2 of ^14j. Let us introduce the weight 



w 



(,,,,2)(x) = (l + x)"^l-a;)"^ -l<x<l (C.28) 



together with the set of polynomials {pf^'"^\x)} orthogonal with respect to ^(0-1,0-2); i-e- 



^K,<x2)(^)pf'^'Hx)p^"^^)(x) da; = S,,. (C.29) 
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We have 

r /_ _ N ' )j,k=0,...,N-l, 

where 



J^(ai, 0^2) = det{a%"'''^),,k=o,...,N-i, (C.30) 



The orthogonal polynomials {pj ^' ^ (x)} are called Jacobi polynomials. 
1. 0"'"(2A^). The average to compute is 



i=l / o+(2Ar) Vi=l 



(C.32) 



^ -P(^l/2 -1/2)(^1' • • • ' ^A^) d6'i ■ ■ ■ d6'Ar. 



The substitution ()(I26|1 gives 



Vi=l / o+(2Af) 

(For simplicity, from now on we shall drop the superscript (ai, O'l) when denoting the 
orthogonal polynomials, their weight and the matrix elements flC.31|) .) When the 
weight is = (1 — a;^)~^/^ the Jacobi polynomials are also known as Chebyshev 
polynomials of the first kind; these are 

1 f2 
Pq{x) = —j= and Pj{x) = y — cos (j cos ^ xj , j > 0- (C.34) 

Substituting x = cos 9 in equation ()C.31|) we have 

1 f^"" 

"00 = 7^ / 9{0) d9 = go, (C.35a) 

aoj = ajo = — [ 9(9) cos(j 9) d9 = V2gj, j > 0, (C.35b) 

Jo 

2 n 

<^jk = - 9{d) cos(j 9) cos(/c 9) d9 = gj^k + 9j+k, j,k>0 (C.35c) 



7r 

1 r27r 



where g, = ^J^'''gi9)e--^^'d9. 
2. Sp(2A^). The average to compute is 



^J=l / Sp(2JV) Vi=l 

^ -P(l/2,1/2)(^1) • • • 5 ^Af) d6'i ■ ■ ■ d6'Ar. 



(C.36) 
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The substitution ()C.26|) gives 



N 



n /(^^•) ) = ^^(V2, 1/2) = det ,=o,...,7v-i • 

^■?'=l ' Sp{2Af) 



(C.37) 



When the weight is w{x) = (1 — x^)^/^ the Jacobi polynomials reduce to the Cheby- 
shev polynomials of the second kind: 



, 2 sin [(? + 1) cos x] 



TT sin fcos ^ X 



(C.38) 



The matrix elements ()C.31|) becomes 

c^jk = - [ 9(0) sin [(j + 1)^] sin [{k + 1)6] dO 

Jo 

= 9j-k- 9j+k+2, j,k = 0, . . . ,N - 1. 



(C.39) 



3. 0^(2A^ + !)• Since the determinant is positive and the dimension is odd, the eigen- 
values come in complex conjugate pairs and the extra eigenvalue is 1. This property 
must be taken into account when computing the average ()(7.24jl . We have the formula 



N 



^•'=1 ' 0+{2Af+l) 



TV •/ —TV 



N 



\j=i 



^ -P{^i/2,i/2)(^i5 • • • 5 ^n) d9i - ■ ■ d^Af, 



which can be rearranged to give 



N 



n /(^^O ) = /(0)/7v(-l/2, 1/2) = /(O) det (a,fc),. ,=o,...,7v- 

^J=l ' 0+(2Af+l) 



1 • 



The Jacobi polynomials with weight w{x) = — x)/{l + x) are 

1 sin[(^)cos-i( 

Pji^) = ri _ 

y/TC sm [2 cos 

The matrix elements ()C.31|) become 



[x 



J>0. 



TT 



^jk = - 5' (6*) sin 



2j + l 



sm 



2k +1 



de 



9j-k-9j+k+i, j,k = 0,...,N -1. 



(C.40) 



(C.41) 



(C.42) 



(C.43) 
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^(<Tl,a2)UV j 








U(iV) 


det {oLjk).^^^^^ 


A^— 1 


/i-fc 


, J,k>0 








90 iij = k = 




det (aifc)j. fc=o,. 


.,N-\ 




\/2gi if either j = 0, k = I 
or j = I, k = 


Sp(2Ar) 


det (ttj-fc)^. fc=o,. 


.,Af-l 


fi'j-fc 


k + gj+k if j,k>0 
- 9j+k+2, j,k>0 


0+(2iV+ 1) 


/(O) det 


0,...,Af-l 




- 9j+k+i, j,k>0 


0-(2iV+ 1) 


/(7r)det 


=0,...,Af-l 




+ gj+k+i, j,k>0 


0-(2iV + 2) 


/(0)/(7r) det 


fc=0,...,Af-l 


9j-k 


- 9j+k+2, j,k>0 



Table 1: Summary of the averages over the different compact groups of F{U) = Y[f=i fi^j)- 
The function / is even and g{9) = f{9)f{—9); fi and gi are the corresponding Fourier coef- 
ficients; is the dimension of the matrices and the number of independent eigenvalues. 



4. O (2A^ + 1). The treatment in this case is similar to that one for 0+(2A^+ 1), except 
that now the extra eigenvalue is —1. Therefore, we have 



N 



n /(^^O ) = /(^)^^( 1/2, -1/2) = /(vr) det {a.^), 

^•?=l ' 0-(2Ar+l) 



j,k=0,...,N-l 



(C.44) 



The weight is w{x) = {1 + x) / {1 — x) and the corresponding Jacobi polynomials 
are 

1 cos [(^) cos-i(x)l , 



y/TT COS [| cos~^(a;j 
The matrix elements ()C.31|) are given by the integral transform 



Oijk 



IT 



g{9) cos 



2j + l 



9 



cos 



2A; + 1 



d^ 



(C.46) 



9j-k + 9j+k+i, j,k = 0,. . .,N -1. 



0~{2N + 2). Since the determinant is negative, one eigenvalue must be —1, and 
since they come in complex conjugate pairs 1 is also an eigenvalue. The number of 
independent eigenvalues is therefore only N and not + 1. The Haar measure is 
given by formula ()C.25|) with {ai,a2) = (1/2, 1/2), the same as Sp(2iV). Therefore, 
we have 



hm)) =/(o)/w/ ■■■/ n/(^^)/M^- 

\J=1 I 0-(2iV+2) Vj=l 

^ -P(l/2,1/2)(^1' • • • ' d6'i ■ ■ ■ d^AT, 



(C.47) 
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which becomes 

JV 



IlM) =/(0)/(vr)/^(l/2,l/2) 

1 / (C.48) 

= /(0)/(vr) det (ajfc)^- fc=o,...,7V-i ' 

where the ajk are given by equation ()C.39|) . 

There exist useful relations among averages of functions over the classical compact 
groups. Those that are relevant to this paper are the following: 

/2N+1 \ /2N+2 \ / 2N \ 

/ U(2Ar+l) \ ' 0+(2Af+2) ' Sp(2Ar) 

HM)) =(n/(^-%)) > (C.49b) 

^■?=1 ' 0+(2Ar+l) ' 0-(2Ar+l) 

flM)) =(n/(^-%)) , (C.49c) 

^■?'=1 ' 0-(2Ar+l) ' 0+(2Af+l) 

2Af \ ^ /2N+2 \ 

n^C'V =7M7w(nw,)) . (c.49d) 
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